Quantum transport in randomly diluted quantum percolation clusters in two 

dimensions 
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We study the hopping transport of a quantum particle through finite, randomly diluted percolation 
clusters in two dimensions. We investigate how the transmission coefficient T behaves as a function 
of the energy E of the particle, the occupation concentration p of the disordered cluster, the size 
of the underlying lattice, and the type of connection chosen between the cluster and the input and 
output leads. We investigate both the point-to-point contacts and the busbar type of connection. 
For highly diluted clusters we find the behavior of the transmission to be independent of the type of 
connection. As the amount of dilution is decreased we find sharp variations in transmission. These 
variations are the remnants of the resonances at the ordered, zero-dilution, limit. For particles 
with energies within 0.25 < E < 1.75 (relative to the hopping integral) and with underlying square 
lattices of size 20x20, the configurations begin transmitting near pa = 0.60 with T against p curves 
following a common pattern as the amount of dilution is decreased. Near Pfs — 0.90 this pattern is 
broken and the transmission begins to vary with the energy. In the asymptotic limit of very large 
clusters we find the systems to be totally refiecting except when the amount of dilution is very low 
and when the particle has energy close to a resonance value at the ordered limit or when the particle 
has energy at the middle of the band. 

PACS numbers: 05.60.Gg,05.50.-|-q,05.10.-a,73.23.-b 



I. INTRODUCTION 

In classical systems a particle can traverse around a 
configuration of barriers as long as there is an energeti- 
cally available path. However, in systems where quantum 
mechanical effects cannot be ignored due to the wavelike 
nature of the particle, constructive and destructive inter- 
ference can occur whenever there are paths of different 
lengths available to the particle. The transport of a par- 
ticle is therefore significantly influenced by quantum in- 
terference in systems where quantum mechanical effects 
cannot be ignored. In a previous paper [l| we focused on 
the effects of quantum interference in the hopping trans- 
port of a particle through finite, ordered square lattices. 
In this paper we extend our studies to include the trans- 
port of a particle through finite disordered clusters in 
two dimensions (2D) where the disorder is introduced by 
random dilution. 

The question of whether disordered clusters in 2D are 
always insulating or not is controversial. One-parameter 
scaling theory ^ states that there is no metallic state in 
non-interacting 2D disordered systems. All states should 
be localized in an infinitely large and disordered 2D sys- 
tem at zero temperature. In weak localization, states 
are logarithmically localized, while in strong localiza- 
tion, states are exponentially localized. In the presence 
of disorder, states are either weakly or strongly local- 
ized depending on the amount of disorder. However, re- 
sults from recent experiments on dilute low-disordered 
Si MOSFET and GaAs/AlGaAs heterostructures show 
hints of a metallic behavior and a metal-to- insulator tran- 



sition in 2D disordered systems [1,11,1^. Several sub- 
sequent experiments on different 2D disordered systems 
also produce results hinting at a metal-to-insulator tran- 
sition. For a review of these experiments see Abrahams et 
al. IQ] Furthermore, experiments done on self-assembled 
quantum dots on Ga[Al]As heterostructures 0] also sug- 
gest a metal-to-insulator transition. Since one-parameter 
scaling theory deals only with non- interacting particles, 
it is possible that a metallic behavior can occur as a re- 
sult of the interplay between inter-particle interactions 
and the disorder. There is also the question of whether 
the type of disorder chosen is relevant. One of the classic 
and intensively investigated model of hopping transport 
of non-interacting particles in disordered systems is the 
Anderson model of localization 8]. The particle in this 
model is governed by the tight-binding Hamiltonian 

i (ij) 

where the \i) and |j) are tight-binding basis functions 
centered on sites i and j, respectively, Vij = 1 if i and j 
are nearest-neighbor sites and Vij = otherwise, and the 
on-site energies e^ are randomly chosen in some range 
Cil < W. A particle traverses the lattice by hopping 
from a site to a nearest neighbor site. As the range W is 
increased a transition from conducting to localized states 
occurs for systems of dimension d > 3. However, all 
states are localized for systems of dimension two or below. 
The type of disorder in this model therefore can not give 
rise to the observed metallic behavior in 2D disordered 
systems. 



The disorder in Eq. |^ lies in the on-site energies e^ 
while the underlying lattice is ordered. A variant of the 
Anderson model is the quantum percolation model |flLllfl| 
wherein the on-site energies e^ are held constant while the 
underl ying lattice is a disordered cluster from percolation 
theory [l^] ■ Since the underlying cluster is disordered not 
every nearest neighbor site is available for the particle 
to hop onto. In matrix representation the disorder in 
the Anderson model is located along the diagonal of the 
Hamiltonian while in quantum percolation the disorder 
is manifest at off-diagonal locations. 

There is a long-standing question whether the Ander- 
son model and the quantum percolation model yield the 
same type of localization behavior and thus, belong to the 
same universality class. For disordered clusters in three 
dimensions it is widely agreed that there is a transition 
from extended to localized states in quantum percola- 
tion milllll El d. For disordered clusters in 2D, 
however, there is no clear consensus whether such a tran- 
sition exists in quantum percolation. Several groups of 
researchers including those using the dlog Fade approx- 
imation method |l7l |. real space renormalization method 
|l8| , or studying the inverse participation ratio [l9| found 
a transition from exponentially localized states to non- 
exponentially localized states for the site concentration 
that ranges within 0.73 < Pq < 0.87. There are, how- 
ever, studies that do not find evidence of such a transi- 
tion. Studies by numerically calculating the conductance 
and making use of the one-parameter scaling hypothe- 
sis |2fl], by using the vibration-diffusion analogy _21j, by 
finite-size scaling analysis and by transfer matrix meth- 
ods [22|, and vector recursion technique [23 found no 
evidence of a transition. A study by Inui et al. |2J| 
found all states to be localized except for those with par- 
ticle energies at the middle of the band and when the 
underlying lattice is bipartite, such as a square lattice. 
It is therefore not clear whether the Anderson model and 
the quantum percolation model produce the same type of 
localization behavior and if not, whether quantum per- 
colation may be relevant to the experimentally observed 
conducting behavior in disordered systems in 2D. 

There are three adjustable parameters in quantum per- 
colation: the energy E of the incident particle, the con- 
centration p for the random dilution of the lattice, and 
the size of the system. In this paper we present our stud- 
ies and results by varying the values of all three param- 
eters. 



II. COUPLING THE LEADS TO THE CLUSTER 
AND DETERMINING THE CONDUCTANCE 

In a previous paper [J] we studied the effects of quan- 
tum interference in the transmission of a particle through 
ordered, i.e., undiluted, square lattices. To study the 
transport properties we connected semi-infinite chains to 
the square lattices. The incident particle was then set to 
propagate from one chain, pass through the square lat- 



tice, and then transmit through to the other chain. We 
used the quantum percolation model and set the constant 
on-site energy e^ = in Eq. (Q . Transmission and reflec- 
tion coefficients were then determined from the eigen- 
states of the whole system including the chains. From 
the transmission and reflection coefficients the conduc- 
tance may be calculated by making use of the Landauer- 
Biittiker formalism 25]. For transport through ordered 
square lattices we found transmission and reflection res- 
onances whenever the energy of the incident particle was 
close to a doubly-degenerate eigenvalue of the square lat- 
tice. We also found the type of connection chosen be- 
tween the chains and the square lattice to strongly in- 
fluence the transport characteristics. Here we follow the 
same method described above in determining the trans- 
port properties of disordered clusters. 
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FIG. 1: An illustration of a disordered cluster on a 3x3 lattice 
attached to the input and output leads through point-to-point 
contacts. Beside each site is its unique label. 
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disordered cluster 

FIG. 2: An illustration of a disordered cluster on a 3x3 lat- 
tice attached to the input and output leads through a busbar 
connection. Beside each site is its unique label. Notice that 
in contrast to the point-to-point contacts there are multiple 
connections between the cluster and the leads in a busbar. 

In this paper the disordered clusters are constructed 
using site percolation on a square lattice [ll|. There are 
several ways of attaching the semi-inflnite chains to the 
disordered clusters. Shown in Figs. ^ and |21 are two pos- 
sible ways of attaching the chains. We label the chain 
where the particle is incident from as the input lead and 
the chain where the particle transmit through as the out- 
put lead. A point-to-point contacts type of connection 
is shown in Fig. ^ The input lead is attached to the 
top-leftmost site while the output lead is attached to the 
bottom-rightmost site in the disordered cluster. The par- 
ticle is incident through the input lead, passes through 



the cluster, and then transmits through the output lead. 
If the cluster is not spanning it is not possible to com- 
plete the connection from the input to the output lead. 
The incident particle is then hindered from propagating 
through the cluster. The minimum necessary require- 
ment therefore for the particle to reach the output lead 
is for the cluster to be spanning. In site percolation on a 
square lattice the critical concentration of occupied sites 
is Pc = 0.5927. In quantum percolation, if there is indeed 
a transition from localized to extended states as the con- 
centration of sites is increased then this transition should 
occur above the classical percolation threshold pc- 

A busbar type of connection between the cluster and 
the leads is shown in Fig. [3 The input lead is attached 
to all the sites at the left side of the cluster while the 
output lead is attached to all the sites at the right side 
of the cluster. In contrast to the point-to-point contacts, 
there are multiple connections between the leads and the 
cluster in a busbar. A physical realization of a busbar 
can be a set-up wherein the contacts connecting the dis- 
ordered cluster to the current source are large enough to 
encompass and attach to the whole side of the cluster. 

There are other ways of connecting the disordered clus- 
ter to the input and output leads. We have, however, cho- 
sen the point-to-point contacts and the busbar because 
at the ordered limit, when the cluster becomes a fully oc- 
cupied lattice, these two types of connections are compli- 
mentary in the sense that point-to-point contacts main- 
tain the bipartite symmetry of the square lattice while 
the busbar, because of its multiple connections, breaks 
that symmetry [l|. 

To determine the transport properties of the disordered 
clusters we follow the method described by Daboul et 



al. [13 Setting the constant on-site energy e^ 
quantum percolation Hamiltonian becomes 
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We then attach the semi-infinite chains to the cluster. 
The transmission and reflection coefhcients can then 
be determined from the resulting eigenvalue equation 
HqpTp = Eip. However, this is an infinitely-sized prob- 



lem and to reduce it to a finite one Daboul et al. 
the following ansatz: 
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where n = 0, 1, 2, . . .. Notice that shown beside each site 
in Figs, n and 121 is its unique label. Sites belonging to 
the cluster are labeled alphabetically. Sites belonging to 
the input chain are labeled by negative integers while 
those belonging to the output chain are labeled by pos- 
itive integers. In the ansatz above ip-(^n+i) is the part 
of the wavcfunction along the input chain and ip^(^n+i) 
is the part of the wavefunction along the output chain. 
The t and r are the transmission and reflection ampli- 
tudes, respectively. What the ansatz implies, therefore, 
is that the incoming plane wave is partially reflected back 
through the input chain and partially transmitted out to 
the output chain. 

Using the ansatz in Eq. |J3J| we can reduce the infinitely- 
sized eigenvalue problem into a finite one. For example, 
for the busbar configuration shown in Fig. |2] the eigen- 
value problem in matrix representation reduces to 
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The connection between the cluster and the leads is man- 
ifest in the perimeter of the 9x9 square matrix above. 
The inner 7x7 inner sub-matrix involves the connections 
among sites in the cluster and is therefore independent of 
the type of connection between the cluster and the leads. 
In order for the ansatz to lead to valid solutions, the fol- 
lowing relation between q and E must be satisfied along 
the chains: 
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The reflection coefBcient i? = |r| and transmission coef- 
ficient T = |i| can then be calculated from Eqs. Q and 
(O and by choosing a value for the incident particle's en- 
ergy E. Note that because of the constraint condition, 
i.e., Eq. (Q, along the chains the particle's energy is now 
restricted to be within —2 < E < 2. Also, although the 
problem now involves a finite matrix, effectively it is still 
an infinite system. Therefore, the energy E is contin- 
uous and any value within the range [—2,2] leads to a 



valid solution. 

Eq. Q is in the form of a linear equation Ax = b with x 
as the unknown. We numerically determine x exactly by 
solving the inverse, A~^, using singular value decomposi- 
tion 26] and then multiplying it to b. For a given lattice 
size, LxL, as we increase the occupation concentration 
p the size of the disordered cluster, on average, increases 
as well. Consequently, the size of the linear problem to 
be solved increases. For a disordered cluster of size TV 
the size of the matrix A to be solved is {N + 2) .In this 
paper the largest clusters we investigate are those embed- 
ded in lattices of size 30x30. We repeat the calculations 
of the transmission amplitude t and reflection amplitude 
r for 1000 realizations of the disorder and then take the 
average of |i| and \r\ to get the average transmission 
coefficient T and reflection coefficient R, respectively. 



III. NUMERICAL RESULTS 

We investigate the characteristics of the transmission 
coefficient by varying the energy E of the incident parti- 
cle, the site occupation concentration p of the disordered 
clusters, and the size LxL of the lattice. Shown in Fig.O 
is a plot of the transmission coefficient T as a function 
of the energy E of the incident particle. The clusters are 
connected to the input and output leads through point- 
to-point contacts. In every plot shown in this paper each 
data point is an average over 1000 realizations of disorder 
configurations. 



cupation concentration goes from p = 0.90 to p = 0.95. 
In particular, at the middle of the band, the transmis- 
sion goes from a dip to a peak. We also start seeing 
more sharp variations in the transmission curve begin- 
ning at p = 0.95. For very low amount of disorder, i.e., 
for p — 0.99, the fluctuations are more pronounced and 
the peak at the middle of the band is near the full trans- 
mission value of T = 1 . In our previous paper |l| we have 
shown that in the case of no disorder, i.e., for p = 1, reso- 
nance occurs whenever the energy of the incident particle 
is close to a doubly degenerate eigenvalue of the square 
lattice without the leads. Thus, the pronounced varia- 
tions that we see in Fig. for low amounts of disorder 
are actually remnants of the resonances that exist at the 
ordered limit. 
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FIG. 3: Plot of the transmission coefficient T as a function 
of the incident particle's energy E for disordered clusters on 
20x20 lattices with point-to-point contacts between the leads 
and the cluster. The disorder concentrations are p — 0.60 
(T), p = 0.70 (V), p = 0.80 (♦), p = 0.90 (0), P = 0.95 (a), 
and p = 0.99 (A). 

In Fig. 13 notice that for p = 0.60, 0.70, and 0.80 when 
the clusters are highly disordered the transmission curves 
are shaped like an inverted w with a dip at the middle 
of the band. The curves then change shape as the oc- 



FIG. 4: Plot of the transmission coefficient T against the 
incident particle's energy E for disordered clusters on 20x20 
lattices with the busbar connection between the leads and 
the cluster. The disorder concentrations are p = 0.60 (T), 
p = 0.70 (V), p = 0.80 (♦), p = 0.90 (0), P = 0.95 (a), and 
p = 0.99 (a). Notice that the transmission values are very 
small. 

Quantum interference occurs within the disordered 
cluster and also at the connections between the leads and 
the cluster. With the busbar, the multiple connections 
between the cluster and the leads enhance the effect of 
quantum interference. An outcome of this is a signifi- 
cant decrease in the transmission of the particle. Shown 
in Fig. 21 is a plot of the transmission coefficient T as a 
function of the energy E of the incident particle for the 
busbar connection. Notice that the scale of the trans- 
mission only goes up to T = 0.06. For highly disordered 
clusters, just like in the point-to-point contacts case, the 
shape of the transmission curve resembles an inverted w 
with the dip at the middle of the band. As the amount 
of disorder is decreased we see a growing asymmetry be- 
tween the positive and negative sides of E. For very 
low amount of disorder, i.e., for p = 0.99, there are pro- 
nounced variations in the transmission. These variations 
are remnants of the resonances at the ordered limit of a 
square lattice connected to leads through a busbar yj]. 



We now investigate how the transmission behaves as 
the disorder concentration is varied while the energy of 
the particle is held fixed. Shown in Fig. is a plot of 
the transmission coefficient T as a function of the site 
occupation concentration p for the point-to-point con- 
tacts type of connection. The results are from the 20 x 20 
lattices. 




FIG. 5: Plot of the transmission coefficient T as a function of 
the cluster occupation concentration p for disordered clusters 
on 20x20 lattices with point-to-point contacts between the 
leads and the cluster. The incident particle has energies E = 
0.00 (♦), E = 0.25 (A), E = 0.50 (a), E = 0.75 (V), and 
E^ 1.00 (T). 

The plot in Fig. [S] is for five selected values of the par- 
ticle's energy E. These values are all positives since as 
we have seen in Fig. |31 there is symmetry between the 
positive and negative sides of E when the connection is 
point-to-point contacts. In Fig. |31 notice that the E — 
case is special. For high disorder concentrations p = 0.50 
until around p = 0.90 the curves for the other E values 
follow the same pattern. The system begins to be trans- 
mitting slightly above p — 0.6. Near p — 0.90 the trans- 
mission curves begin to spread out and at the region near 
p = 1 there are sharp variations as the system seems to 
change abruptly from one dominated by disorder to one 
that reflects the ordered limit wherein resonances occur. 
For the system with E = it does not begin to trans- 
mit until around p = 0.80. The transmission curve then 
goes up with a slope that is steeper than the curves for 
the other E values. In Fig. O we have seen that at the 
middle of the band the transmission shifts from being a 
dip to a peak as p is varied. This shift in the behavior of 
the transmission is also manifest in the rapid rise of the 
transmission curve after p = 0.90 for the E = case in 
Fig.El 

For the busbar connection, we have seen in Fig.Qlthat 
there is no symmetry between the positive and negative 
sides of the energy E. We therefore show in Figs. and [7| 
plots of the transmission coefficient T as a function of the 
occupation concentration p as the energy E of the inci- 
dent particle is fixed at either negative or positive values, 



respectively. The system with the busbars is highly re- 
fiecting and we notice therefore in Figs. El and [7| that the 
transmission scales are only up to T = 0.02. Just like in 
the case for the point-to-point contacts, the E = case 
appears to be special. 
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FIG. 6: Plot of the transmission coefficient T" as a function of 
the cluster occupation concentration p for disordered clusters 
on 20x20 lattices with busbar connections between the leads 
and the cluster. The incident particle has energies E = 0.00 
(♦), E = -0.25 (A), E = -0.50 (a), E = -0.75 (V), and 
£ = -1.00 (T). 
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FIG. 7: Plot of the transmission coefficient T as a function of 
the cluster occupation concentration p for disordered clusters 
on 20x20 lattices with busbar connections between the leads 
and the cluster. The incident particle has energies E = 0.00 
(♦), E = 0.25 (A), E = 0.50 (a), E = 0.75 (V), and E = 1.00 

In Fig. El the selected energies are negative while in 
Fig. [3 the energies are positive. The data points for 
E — are shown in both figures. Notice again that 
these points do not follow the same pattern as the data 
points for the other energies. In particular, the system 
is non-transmitting for high disorder concentrations until 



near p = 0.90. For the other E values the system begins 
to transmit near p — 0.60. After about p = 0.90 the 
transmission curves spread out just hke in the case for 
point-to-point contacts. 

For both the busbar and point-to-point contacts we 
thus see signatures of a shift in transmission characteris- 
tics near pa = 0.60 and then again near pp = 0.90. The 
system begins to transmit near pa regardless of the type 
of connection chosen between the cluster and the leads, 
except at the middle of the band. Note that this value of 
p is close to the classical percolation threshold Pc of site 
percolation on a square lattice. In order for the parti- 
cle to transmit through a disordered cluster that cluster 
must at least be spanning. In an infinite system, a span- 
ning cluster only appears above the classical Pc- How- 
ever, due to quantum interference and the finite sizes of 
the systems we study, the configuration will not neces- 
sarily begin transmitting exactly at pc- The shift near pp 
signals the onset of the dependence of the transmission 
coefficient T on the particle's energy E as the occupation 
concentration p is held fixed. Between pa and pp the 
transmission behaves in the same fashion, except, again, 
at the middle of the band. This can also be seen in Fig.O 
where for a high disorder concentration, say for p = 0.80, 
the transmission along the interval 0.25 < E < 1.75 is 
relatively flat. 

We further investigate the characteristics of the trans- 
mission as the size of the lattices are varied. Shown in 
Figs. IHl and |51 are plots of the transmission as a function 
of lattice sizes as the occupation concentration p and en- 
ergy E of the incident particle are both held fixed. Also 
shown in the plots are the best-fitting exponential curves 
and best-fitting power-law curves. 





FIG. 9: The transmission coefficient T as a function of the 
lattice length L with an incident particle of energy E — 0.00. 
The occupation concentrations are p — 0.80 (T), p = 0.92 
(a) and p — 0.99 (♦). The solid lines are the best-fitting 
exponential curves while the dashed lines are the best-fitting 
power-law curves. For p = 0.80 the three points after L = 25 
have transmissions that are practically zeroes and are there- 
fore not included when the best-fit lines are determined. 



'AS E — 1. We have also chosen the connection between 
the cluster and the leads to be point-to-point contacts be- 
cause this connection type is generally more transmitting 
than the busbar. The cluster occupation concentrations 
are chosen as p = 0.68 and p = 0.92 since these values 
are near pa and pp. The power-law fits are in the form 
T — TpL^'^p while the exponential fits are in the form 
T = Tee-'^"^, for lattices of size LxL. Shown in Tabled 
are the fitting values corresponding to the best-fit curves 
in Fig. |S1 The Rp and Re are the correlation coefficients 
for the power-law and exponential fits performed as linear 
regressions of log (T) against log (L) and log (T) against 
L, respectively. 
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TABLE I: Table of fitting values for the E ^ 1 case. The Tp 



(jp, and Rp belong to the power-law fits while the Te 
Re are for the exponential fits. 
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FIG. 8: Plot of the transmission coefficient T as a function 
of the lattice length L. The energy of the incident particle is 
E = 1.00. The cluster and the leads are connected through 
point-to-point contacts. The occupation concentrations are 
p = 0.68 (A), p = 0.80 (♦), and p = 0.92 (t). The solid lines 
are the best-fitting exponential curves while the dashed lines 
are the best-fitting power-law curves. 

In Fig. |S1 we have chosen the energy of the particle 



From the values of the correlation coefficients, Rp and 
Re, it is not possible to clearly distinguish between the 
goodness of the power-law and exponential fits as they 
are very close to each other. However, since these cor- 
relation coefficients are reasonably close to 1, we can be 
fairly confident that T extrapolates to zero as L — > oo 
or in the thermodynamic limit. Thus, in either case, our 
result for ii^ = 1 is consistent with one-parameter scal- 
ing theory |3| which predicts that wavefunctions in two 



dimensions are localized. 

Shown in Fig. |51is a plot of T as a function of L when 
the particle's energy \s E = Q. The corresponding best-fit 
estimates are shown in Tablellll While the power-law and 
exponential curves are both reasonably good for p = 0.80, 
the power-law fit is not nearly as reasonable for p = 0.92. 
For very low amount of disorder at p = 0.99 neither the 
power-law nor the exponential fit is reasonable. From 
the values of a^p and CToe we get a hint that the T can 
be independent of L at very low disorder concentrations, 
at least at certain values of E. For such a situation, it 
is possible for the transmission to be non-zero even for 
asymptotically large clusters. 
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TABLE II: Table of fitting values for the £ = case. The 
Tbp, o"op, and Rop belong to the power-law fits while the Toe, 
(Toe, and Roe are for the exponential fits. Note that neither 
fits are good for the p = 0.99 case. 



IV. SUMMARY AND CONCLUSIONS 

The input and output leads can be attached to the 
disordered cluster in several different ways. In this paper 
we have chosen the connections to be point-to-point con- 
tacts and the busbar. The behavior of the transmission 
T is independent of the type of connection chosen for 
highly disordered clusters. This is because the number 
of actual multiple connections in the busbar case is not 
large enough to differentiate its effects from that of the 
point-to-point contacts type of connection. As the site 
dilution is diminished, i.e., as p increases, the effects of 
the choice of connection begin to be apparent. For low 
disorder, i.e., high p, point-to-point contacts, in general, 
are more transmitting than the busbar. The resonances 
at the ordered limit, however, influence the transmission 
at low disorder. We thus see transmission variations near 
p = 1 as we vary the energy of the particle. 

We have shown how the transmission coefficient T be- 
haves as a function of the energy E of the particle, the 
occupation concentration p of the disordered cluster, and 
the size LxL oi the underlying lattice. For 20x20 lat- 
tices we find the system to begin transmitting around 
Pa = 0.60, which is near the classical percolation thresh- 
old pc- The transmission is independent of the incident 



particle's energy along the interval 0.25 < E < 1.75 when 
p is increased until around Pf3 — 0.90. After pp the trans- 
mission begins to be dependent on the energy of the par- 
ticle. The exception to this pattern is when the particle 
has energy at the middle of the band. For such a particle 
the system does not begin to be transmitting until after 
around p = 0.80 in point-to-point contacts and p = 0.90 
in the busbar. 

Extrapolating the behavior of the transmission to very 
large lattices, we find T — > in moderate to high disor- 
der. The system is localized in these situations. The be- 
havior of the transmission at very low disorder, however, 
is not as straightforward. The resonances at the ordered 
limit influences the behavior of the transmission at very 
low disorder. So even with asymptotically large clusters, 
as long as the disorder is low, it may be possible to get 
non-zero transmission when the energy of the incident 
particle is close to the resonances at the ordered limit. 
When the energy of the incident particle is at the mid- 
dle of the band, the system can also be transmitting at 
very low disorder even for asymptotically large clusters. 
This case is consistent with the work of Inui et al. |2J| 
where they cite the bipartite symmetry of the underly- 
ing square lattice as the reason for finding wavefunctions 
that are not exponentially localized at the middle of the 
band. 

Quantum percolation in two dimensions in asymp- 
totically large clusters therefore generally lead to non- 
conducting systems. Our results suggest that the only 
possible exceptions occur at very low disorder when the 
incident particle's energy is either at the middle of the 
band or near the resonance value in the ordered limit. 
Quantum percolation is a single-particle quantum model 
of hopping transport where the effects from mechanisms 
such as tunneling, long-ranged hopping, or inter-particle 
interactions are not taken into account. Taking these 
other effects into account may enhance the transmission 
of the particle in such a way that the configuration may 
be transmitting even for asymptotically large and highly 
disordered clusters. 
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